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Abstract 

Background: Bacterial infections of wounds impair healing and worsen scarring. We hypothesized that 
transcriptome analysis of wounds infected with Klebsiella pneumoniae (K.p.) or Pseudomonas aeruginosa (P.a.) would 
indicate host-responses associated with the worse healing of P.a- than K.p -infected wounds. 

Methods: Wounds created on post-operative day (POD) 0 were infected during the inflammatory phase of healing 
on POD3 and were harvested on POD4 for microarray and transcriptome analysis. Other wounds received topical 
antibiotic after infection for 24 hours to promote biofilm development, and were harvested on POD6 or POD12. 

Results: Wounds infected for 24 hours, relative to uninfected wounds, elevated transcripts of immune-response 
functions characteristic of infiltrating leukocytes. But P.a.-infected wounds elevated many more transcripts and 
to higher levels than Kp.-infected wounds. Coincidently, suppressed transcripts of both wounds enriched into 
stress-response pathways, including EIF2 signaling; however, this was more extensive for P.a.-infected wounds, 
including many-fold more transcripts enriching in the 'cell death' annotation, suggesting resident cutaneous cell 
toxicity in response to a more damaging P.a. inflammatory milieu. The POD6 wounds were colonized with biofilm 
but expressed magnitudes fewer immune-response transcripts with no stress-response enrichments. However, 
elevated transcripts of P.a.-infected wounds were inferred to be regulated by type I interferons, similar to a network 
unique to P.a.-infected wounds on POD4. On POD12, transcripts that were more elevated in Kp.-infected wounds 
suggested healing, while transcripts more elevated in P.a.-infected wounds indicated inflammation. 

Conclusions: An extensive inflammatory response of wounds was evident from upregulated transcripts 24 hours 
after infection with either bacterium, but the response was more intense for P.a.- than Kp.-infected wounds. 
Coincidently, more extensive down-regulated transcripts of P.a.-infected wounds indicated a stronger "integrated 
stress response" to the inflammatory milieu that tipped more toward cutaneous cell death. Unique to P.a.-infected 
wounds on POD4 and POD6 were networks inferred to be regulated by interferons, which may result from intracellular 
replication of P.a. These data point to specific downregulated transcripts of cells resident to the wound as well as 
upregulated transcripts characteristic of infiltrating leukocytes that could be useful markers of poorly healing wounds 
and indicators of wound-specific treatments for improving outcomes. 
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Background 

Wound healing proceeds in an orchestrated, highly ordered, 
and overlapping cascade of hemostasis, inflammation, prolif- 
eration, and remodeling [1,2] that can be deranged by local 
or systemic pathologies such as stress [3,4], diabetes [5-9], 
and infections [10-13], including biofilm infections [14,15], 
which can result in delayed wound closure and worsened 
scarring. 

Microbial colonization of both acute and chronic 
wounds is inevitable. In many settings, gram-positive en- 
dogenous skin flora predominate [16]. While this is also 
the case for combat wounds [17], combat casualties hos- 
pitalized for definitive tertiary care in medical treatment 
facilities are at high risk for nosocomial infections 
[18,19] that develop days after injury and are largely 
due to multi-drug resistant gram-negative organisms in- 
cluding Acinetobacter, Pseudomonas, Enterobacter, and 
Klebsiella [18,20-22]. 

Biofilms are highly differentiated and spatially orga- 
nized three-dimensional structures consisting of matrix- 
enclosed communities of one or more microbial species 
that form on colonizable surfaces [23]. Biofilms form 
rapidly in acute wounds [24,25] and contribute to the 
pathogenesis of chronic non-healing wounds [26-28]. 

Inflammation is the normal response to tissue injury 
and microbial infection, and in proper degree is required 
for tissue repair and healing. However, persistent inflam- 
mation or uncontrolled inflammation such as observed 
in chronic and/or infected wounds results in poor 
wound outcomes due to delayed wound closure and 
hypertrophic scarring [29]. The extent of inflammation 
and its clearance in a timely fashion can influence the 
rate and quality of wound healing [30,31]. 

Differences in the host-response to pathogens have been 
associated with virulence mechanisms. P. aeruginosa is a 
common gram-negative bacterium of nosocomial and life- 
threatening infections of immuno-compromised patients 
[32]. It possesses many virulence factors such as exopro- 
teases (e.g., elastase), siderophores, exotoxins, hydrogen 
cyanide, and pyocyanin to attack host defenses. K. pneu- 
moniae account for a substantial percentage of nosocomial 
infections in neonates, patients undergoing respiratory 
therapies, and patients hospitalized in urology and burn 
wards, but less is known about its virulence in wound in- 
fection and healing. Factors such as MagA (responsible for 
the mucoviscosity that correlates with high serum resist- 
ance in vivo [33]), siderophore aerobactin [34], and cell 
surface-associated fimbriae [35] have been implicated in 
K. pneumoniae virulence. 

Using a full-thickness dermal rabbit ear wound model, 
we previously demonstrated that bacterial wound infec- 
tions transitioned from active planktonic-phase infections 
to biofilm-dominant infections that delayed wound clos- 
ure (epithelialization) and granulation tissue in-growth 



[11-13]. The extent of wound healing impairment de- 
pended on the infecting bacterial species. Pseudomonas 
aeruginosa (P.a.) caused the most healing impairment, 
Klebsiella pneumoniae (Kp) slightly delayed closure rela- 
tive to uninfected wounds, and Staphylococcus aureus was 
intermediate between P.a. and Kp. [12]. These infections 
stimulated the expression of inflammatory cytokines such 
as interleukin l-(3 (IL-1|3) and tumor necrosis factor-a 
(TNF-a) to levels proportional to the wound healing im- 
pairment [12]. 

Global gene expression analysis using microarrays has 
been used to define events in tissue repair [1], such 
as those that differ between healing and non-healing 
venous leg ulcers [36], skin wounds and oral mucosa 
wounds (a "privileged" tissue that heals rapidly without 
scar) [37], as well as fetal wounds which heal without 
scar and postnatal wounds that scar [38] . 

Here, we have used microarray/transcriptome analyses 
of wounds infected with P.a. or Kp. to identify similar- 
ities and differences in the wound responses to these 
differently virulent bacteria. The findings contribute to 
understanding the molecular foundations of the im- 
paired healing caused by bacterial infections. 

Methods 

Animals 

Under an approved protocol by the Animal Care and Use 
Committee at Northwestern University, adult female 
New Zealand white rabbits (3-6 months, ~3 kg) were 
acclimated to standard housing and fed ad libitum. All 
animals were housed in individual cages under constant 
temperature and humidity with a 12-hour light-dark cycle. 

Bacterial strains and culture 

P.a. strain PAOl was obtained from the laboratory of 
Dr. Barbara H. Iglewski, University of Rochester Medical 
Center. Kp. strain BAMC 07-18 was kindly provided by 
LTC Clinton Murray of Brooke Army Medical Center, 
Fort Sam Houston. To prepare the bacteria for the 
wound infection, P.a. PAOl and Kp. BAMC 07-18 were 
grown in Luria (LB) broth and tryptic soy broth (TSB), 
respectively. After overnight incubation at 37°C, 0.5-mL 
of the P.a. and Kp. cultures were transferred to 9.5-mL 
of fresh sterile LB and TSB, respectively, and incubated 
at 37°C until the cultures reached the log phase. Bacteria 
were harvested and washed in phosphate-buffered saline 
(PBS) once by centrifugation at 4,000 rpm for 15-minutes 
at 4°C. The resultant pellets were resuspended in PBS and 
adjusted to an optical density of 0.5 at the wavelength 
of 600-nm (OD 600 ). An OD 600 0.5 was equivalent to 
10 8 CFU per mL, which was determined empirically for 
each strain of bacteria used. Ten-|iL of bacteria suspen- 
sion at OD 60 o 0.5 from each strain was used as the 
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inoculum (approximately 1 x 10 6 CFU per wound) for in- 
fecting the rabbit ear wounds. 

Wounding and infection 

Wounding, bacterial infection, and biofilm formation 
were adopted from our previously published wound bio- 
film model [11-13]. Briefly, rabbits were anesthetized 
with an intramuscular injection of ketamine (22.5 mg/ 
kg) and xylazine (3.5 mg/kg) mixture prior to surgery. 
Ears were shaved, sterilized with 70% ethanol, and 
injected intradermally with 1% lidocaine / 1:100,000 
epinephrine at the wound sites. A total of 18 rabbits 
were used. Four, 6-mm diameter, full-thickness dermal 
wounds were created in the central area of the ventral 
ear down to perichondrium (4 wounds per ear). To 
minimize variations in dermal thickness and blood flow 
to the wounds, distal and proximal areas of the ear 
were not used. After surgery, wounds were dressed with 
Tegaderm (3 M Health Care, St. Paul, MN), a semi- 
occlusive transparent film. Individual wounds were inoc- 
ulated with 10 6 bacteria on postoperative day (POD) 3 
and redressed with Tegaderm. Bacteria were allowed to 
proliferate in vivo under the Tegaderm dressing. A top- 
ical antibiotic, Ciloxan ointment (Ciprofloxacin 0.3%, 
Alcon, Fort Worth, TX), was applied on POD4 to elim- 
inate free-floating, planktonic-phase bacteria, leaving a 
predominate biofilm. To prevent seroma formation and 
re-growth of planktonic bacteria, an antimicrobial, ab- 
sorbent dressing containing polyhexamethylene biguan- 
ide (Telfa AMD, Tyco Healthcare Group, Mansfield, 
MA) was applied to wounds on PODs 5, 6, and then 
every other day until harvest. Uninfected and infected 
wounds were equally treated with the antibiotics and 
dressings. All dressings were examined daily throughout 
the protocol. 

Harvesting of wounds 

After euthanizing the animals by intracardiac euthasol 
injection, wounds were harvested for RNA extraction. 
All wounds were excised using a 7-mm biopsy punch 
(Acuderm inc., Fort Lauderdale, FL). Wounds from dif- 
ferent study groups were harvested at POD 4 (prior to 
the application of local antibiotic), POD6, and POD 12, 
and RNA was isolated for microarray and reverse tran- 
scription quantitative PCR (RT-qPCR) analyses. 

RNA extraction 

For mRNA extraction, the dermal layer on the dorsal 
side of the ear wound was removed and the wound bed 
was punched out and samples of all wounds were frozen 
in RNAlater (Ambion, Austin, TX) at -80°C for molecu- 
lar analysis. Wound samples were homogenized using a 
Mini-bead beater-8 (Biospec Products Inc, Bartlesville, 
OK) using Zirconia beads (2.0 mm diameter, Biospec 



Products Inc) in the presence of Trizol Reagent (Sigma- 
Aldrich, St. Louis, MO). Total RNA was isolated according 
to the manufacturers protocol. Contaminating genomic 
DNA during RNA preparation was removed using the 
Turbo DNA- free kit (Ambion, Austin, TX). In general, 
1 (ig of RNA was obtained from each wound. The quantity 
and quality (Absorbance 260/Absorbance 280 ratio, 28 
s/18 s RNA ratio, RNA integrity number) of the extracted 
RNA was determined using a NanoDrop 2000C UV-vis 
Spectrophotometer (Thermo Scientific, Wilmington, DE) 
and the Agilent 2100 Bioanalyzer (Agilent Technologies, 
Santa Clara, CA), respectively. 

Microarray 

RNA (1 \ig each) was collected from 5 control (unin- 
fected), 5 K.p. -infected, and 5 P. <2. -infected wounds at 
each time point (POD4, 6, and 12; 45 samples total). 
Each RNA sample was labeled with the One- Color 
Quick- Amp Labeling Kit (Agilent Technologies) and hy- 
bridized to an Agilent custom rabbit microarray (Agilent 
Product # G2519F) containing 4 sets of 44,000 probes 
per slide. The microarray contains 43,803 rabbit probes 
representing a total of 19,938 unique rabbit sequences. 
All labeling and hybridization was performed at the Uni- 
versity of Florida Interdisciplinary Center for Biotechnol- 
ogy Research according to standard Agilent recommended 
protocols. 

Data collection and analysis 

Data preprocessing was performed using the Bioconduc- 
tor Limma package. The background correction was per- 
formed using Normexp with offset value 50. Quantile 
normalization was then used to correct for the non- 
biological variations across arrays. Probes were filtered 
based on the criteria that at least 39 samples must have 
a signal level above background and negative controls. 
The signal intensity values were further normalized by 
subtracting (in log2 scale) the time zero sample mean 
(i.e., POD3, immediately after wounding). Differentially 
expressed genes were identified using routines imple- 
mented in the Limma package to fit linear models to the 
normalized expression values [39]. The variance used in 
the t-score calculation was corrected by an empirical 
Bayesian method for better estimation under small sam- 
ple size. All ^-values were corrected by the Benjamini- 
Hochberg procedure. 

Ingenuity pathway analysis 

Ingenuity Pathway Analysis (IPA, Ingenuity Systems Inc.) 
was used to analyze enrichments of DEGs (infected vs. un- 
infected wounds, and K.p.- vs. P. <z. -infected wounds) into 
functions, pathways and networks. Because IPA does not 
analyze rabbit gene identifiers (IDs), they were first con- 
verted to their human ortholog gene IDs. Out of the 



Leung et al. BMC Clinical Pathology 2014, 14:20 
http://www.biomedcentral.eom/1 472-6890/1 4/20 



Page 4 of 14 



19,879 sequences on the microarray, 13,481 unique se- 
quences mapped to human IDs and were 'Analysis-Ready 
in IPA (68%). Data were analyzed through the use of 
causal analysis approaches in IPA to infer upstream bio- 
logical causes and probable downstream biological effects, 
using enrichment score (Fishers Exact Test ^-value) to 
measure overlap between observed and predicted gene 
sets, and z-score to match observed and predicted up/ 
down regulation patterns [40] . 

Reverse-transcription qPCR 

Subsets of RNAs were used for determining transcript 
levels of specific target genes by RT-qPCR. RNA was 
used to prepare cDNA using superscript II (Invitrogen) 
with 100-ng of random primers (Invitrogen). RT-qPCR 
analyses using SYBR green was performed utilizing an 
ABI prism 7000 sequence detection system (Applied 
Biosystems, Foster City, CA). Genes that were highly 
differentially regulated in either K.p. or P.a. infected 
wounds as shown in the microarray analyses were selected 
for validation by RT-qPCR. PCR primers were designed 
using the Primer 3 program (http://frodo.wi.mit.edu/). Ex- 
pression of each gene was normalized to the level of the 
house keeping gene, glyceraldehyde 3-phosphate dehydro- 
genase (GAPDH), to get ACt. The 2" AACt method was used 
to calculate gene expression. 

Results 

To identify rabbit ear wound transcriptomic responses 
to infections with K. pneumoniae {K.p.) or P. aeruginosa 
{P.a.)> we used the experimental design summarized in 
Figure 1. Six-mm diameter wounds were created and in- 
fected 3 days later (post-operative day three, POD3) with 
planktonically grown K.p. or P. a. (10 6 colony- forming 
units, CFU). At three post-infection time points, wound 
RNA was harvested for microarray analysis. Wounds 
were first harvested the day after infection on POD4. 
Other wounds infected for 24 hours were treated with 
local antibiotic on POD4 to kill off planktonic bacteria. 
The RNA was isolated from these biofilm-infected 



wounds [12] two days later on POD6, and on POD 12, 
when uninfected wounds were fully epithelialized, K.p.- 
infected wounds were slightly open (-1.7 mm epithelial 
gap [12]), and P. <2. -infected wounds were largely open 
(~4 mm epithelial gap [12]). Gene expression levels of 
the infected wounds at each time point were compared 
with uninfected wounds, or the two types of infected 
wounds were compared each other, and the differentially 
expressed transcripts (DEGs) were analyzed for associa- 
tions of with pathways, networks, and biological func- 
tions and disorders using the Ingenuity Pathway Analysis 
(IPA) database and software. 

Tallies of differentially expressed genes (DEGs) in infected 
vs. uninfected wounds 

We found 1,378 differentially expressed sequences (fold 
change > 2 with adjusted j?-value < 0.05) in infected as 
compared to uninfected (I/U) wounds out of the 19,879 
unique sequences on the microarray; i.e., 7% of all gene 
probes were I/U differentially expressed (I/U-DEGs) by 
>2-fold, when tallying across all PODs for both P.a.- and 
K.p.- infected wounds (Table 1). 

Averaging over all PODs, P.a.- infected wounds con- 
tained 3-fold more I/U-DEGs than K.p. -infected wounds. 
This was largely due to the 5-fold more I/U-DEGs in 
P.a.~ than K.p. -infected wounds on POD4 the day after 
bacteria were inoculated into the wounds and before 
antibiotic treatment [12]. At this time exudate was abun- 
dant and the number of >2-fold-changed I/U-DEGs was 
219 for Kj?.-infected wounds and 1,113 for P. <2. -infected 
wounds (Table 1). There were many fewer DEGs on 
POD6 {K.p.: 17; P.a.: 30) following two days of antibiotic 
treatment when ~1 to 2 logs fewer CFU were recovered 
than were inoculated (our previous study [12]). On 
POD12, the I/U-DEGs had increased {K.p.: 162; P.a.: 
161, Table 1). 

The complete list of all 19,879 gene probes and associ- 
ated normalized log2 fold-changes, adjusted j?-values, 
and updated annotations is provided in Additional file 1, 
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Figure 1 Study design. On postoperative day 0 (PODO), 6-millimeter diameter wounds were created in the central area of the ventral ear down 
to the perichondrium. On P0D3, wounds were inoculated with either P. oeruginoso or K. pneumoniae (10 6 CFU/wound). On POD 4, topical 
ciprofloxacin ointment was applied. Antimicrobial absorbent dressing changes were made on PODs 5, 6 and every second day until wound 
harvest. Wounds were harvested on PODs 4, 6 and 12 by excision using a 7-millimeter biopsy punch, and RNA was isolated for microarray analysis. 
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Table 1 Summary of genes with >2-fold differential 
expression in K.p.- or P.a.-infected wounds vs. uninfected 
wounds 





POD4 




POD6 




POD12 




Infecting bacteria: 


K.p. 


P.a. 


K.p. 


P.a. 


K.p. 


P.a. 


Upregulated genes*: 


168 A 


701 


12 


21 


77 


129 


Down regulated genes*: 


51 


412 


5 


9 


85 


32 


Totals 


219 


1113 


17 


30 


162 


161 



*Based on > 2 fold change of probe intensity with an adjusted p-value < 0.05. 
""Overlapping l/U-DEGs are included. 



and these data have been deposited in the NCBI GEO 
database (Accession number GSE51167). 

Few l/U-DEGs overlapped between PODs, but many 
l/U-DEGs were common to K.p.- and Ra.-infected 
wounds on each POD 

To compare overlapping genes between PODs and infec- 
tion types on each POD, we used the top 200 differentially 
expressed genes, consisting of the 100 most up-regulated 
(Upr) and 100 down-regulated (Dnr) I/U-DEGs, in order 
to compare equal numbers of I/U-DEGs. Not a single gene 
overlapped between all three PODs, for either species, and 
less than 7% overlapped between any two PODs, indicat- 
ing the largely different states of infection and healing 
of the wounds on the three PODs. In contrast, many 
I/U-DEGs were common to K.p.- and P. ^.-infected 
wounds (Figure 2) indicating similar host-responses to the 
different infections. To go beyond comparing similarities 
and differences of single genes, we analyzed similarities 
and differences of pathways and functions associated with 
the I/U-DEGs common and unique to the two infected 
wound types. 



POD4: Upr-l/U-DEGs common to K.p.- and Ra.-infected 
wounds enriched in inflammatory-response functions 

The common I/U-DEGs to K.p.- and P. <z. -infected 
wounds enriched in the IPA Function Annotation cat- 
egory Inflammatory Response'. Downstream Effects 
Analysis of these common genes identified Functions 
Annotations whose activation state was inferred to in- 
crease, based on literature compiled in the Ingenuity® 
Knowledge Base, to include: adhesion of immune cells', 
synthesis of reactive oxygen species', recruitment of neu- 
trophils', activation of myeloid cells', activation of leuko- 
cytes', production of reactive oxygen species', adhesion of 
mononuclear leukocytes; activation of phagocytes', mi- 
gration of neutrophils', cell movement of myeloid cells', 
and 'inflammatory response'. The Z-scores for enrich- 
ment of the I/U-DEGs in these functions were between 
3.2 and 2.4 (significance cut-off 2.0), predicting an acti- 
vated state. The top four inferred upstream regulators for 
these Upr I/U-DEGs common to both infected wound 
types were lipopolysaccharide (LPS), TNF, IFNG and IL1B 
(z-scores 3.5 to 2.5). These results show that the Upr I/U- 
DEGs common to the two infected wound types were 
mostly involved in the inflammatory response. 

POD4: 99% of I/U-DEGs were concordantly Upr or Dnr in 
K.p.- and Ra.-infected wounds 

Seeking to identify differences between the host-responses 
to the two infected wound types, we mined more deeply 
into the dataset, analyzing all I/U-DEGs on POD4 with a 
p- value of <0.05 (Benjamini-Hochberg). The I/U-DEGs 
common to both K.p.- and P.a. -infected wounds summed 
to 907 (Le., all overlapping I/U-DEGs). All of these I/U- 
DEGs, except for 1%, or 9, showed the same direction of 
differential expression in both infected wounds relative to 




0% -\— — i— 

4 6 12 

Post-Operative Day 



Figure 2 I/U-DEGs in common between K.p.- and P.a.-infected wounds on POD4. Percent of Upr (red) and Dnr (green) I/U-DEGs in common 
between wounds infected by P.a. and K.p. on POD 4, 6 and 12 is plotted. The 100 most Upr and 100 most Dnr genes were analyzed for overlaps 
at each POD. Numbers in parentheses are the common/total I/U-DEGs between the two infected wound types. Fewer than 100 Upr and 100 Dnr 
I/U-DEGs were present on POD6. 
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uninfected wounds (45% Upr and 54% Dnr). This result 
suggests biological significance of even small I/U fold- 
changes and indicates largely similar responses of the 
wounds to the two types of bacterial infection. 

The great majority of Upr-transcripts, 91%, were up- 
regulated to a greater extent in P.^.-than K.p. -infected 
wounds. Many of these are mainly expressed in leuko- 
cytes, for example: the superoxide-producing NADPH 
oxidase expressed in neutrophils, NCF1, 2 and 4 (2.8-, 
3.8- and 2.9-fold); the G-protein coupled receptor EMR2 
(2.9-fold); S100A9 (1.9-fold); CSF3R (colony stimulating 
factor 3 receptor, granulocyte (2.9-fold); CSF2, colony 
stimulating factor 2, granulocyte-macrophage (5-fold); and 
CSF3, colony stimulating factor 3, granulocyte (3.6-fold). 
Conceptually, as compared to K.p. -infected wounds, some 
elevated transcripts of P. <2. -infected wounds could come 
from greater numbers of infiltrating leukocytes. Similar 
numbers of leukocytes or cutaneous cells could have up- 
regulated transcription to a greater extent in the absence of 
greater immune cell infiltration, but this seems less likely. 

The majority of Dnr I/U-DEGs, 73%, were also down- 
regulated to a greater extent in P. <z. -infected wounds. 
More Dnr genes and a greater extent of downregulation 
of most overlapping I/U-DEGs in P. a. -infected wounds 
could represent a stronger stress response of cutaneous 
cells to greater inflammation caused by P. a. infection. 

POD4: Transcripts upregulated to a greater extent in 
/C.p.-infected wounds appear to come from cutaneous 
cells responding to the inflammatory milieu 

Of the 9%, or 36, of the concordantly Upr-transcripts 
that were upregulated to a greater extent in /^.-infected 
wounds, only KRT8, CYR61 and CD5 were enriched in 
the function 'inflammatory response'. However, KRT8 
and CYR61 are likely not from leukocytes. KRT8 (kera- 
tin 8) is a type II intermediate filament that forms 
intermediate-sized filaments in the cytoplasm of epithe- 
lial cells [39]. And in fibroblasts, CYR61 appears to up- 
regulate the expression of a number of genes involved in 
angiogenesis, inflammation, and matrix remodeling dur- 
ing wound healing [39]. The top 'Physiological System 
Development and Function' was 'Connective Tissue Devel- 
opment and Function' with seven enriching genes (CYR61, 
IQGAP1, NF2, KRT8, ALDH1A3, MXD1, and CYP27B1). 
Also another keratin, KRT78, was among the 36 genes. 
These enrichments suggest the genes more upregulated in 
K.p. -infected wounds came from cutaneous cells. 

The second topmost 'Molecular and Cellular Function' 
into which the 36 Upr genes enriched was 'Cell Death 
and Survival'. Eleven genes enriched in this function 
(CD5, AIFM3, CYR61, EGLN3, ALDH1A3, KRT8, 
MTFP1, MXD1, NF2, PDK1, and SLC5A8). The topmost 
Upr genes to a greater extent in K.p.- than P.a. -infected 
wounds were CYR61, EGLN3, and HILPDA. These three 



genes are all hypoxia inducible. Additionally, the top ca- 
nonical pathway enrichment of the 36 genes was glycolysis. 

The 36 genes that are Upr to a greater extent in K.p.- 
infected wounds are listed in Additional file 2. Their lack 
of extensive enrichment in inflammation functions and 
their enrichment in functions related to surviving the 
inflammatory environment suggest they represent the 
cutaneous cell response to the inflammatory milieu. Up- 
stream regulators with a predicted activated state for 
regulating these 36 genes were only two: HIF1A and 
NUPR1. HIF1A is the master regulator of the response 
to hypoxia; under hypoxia it activates transcription of 
over 40 genes [41]. NUPR1 is as a component of the 
NUPR1/RELB/IER3 survival pathway that binds to chro- 
matin to change gene expression that is associated with 
resistance to stress [41] and inhibition of apoptosis [42]. 

As an additional analysis, we directly compared differ- 
ential expression between K.p.- and P. <2. -infected wounds 
(i.e., K/P-DEGs), as opposed to the previous analyses 
of differential expression relative to uninfected wounds 
(i.e., I/U-DEGs). The K/P-DEGs that were more elevated 
in K.p. -infected wounds summed to 836 (p-value <0.05, 
Benjamini-Hochberg). The topmost four canonical path- 
way enrichments of these were 'Cell Cycle Control of 
Chromosomal Replication', 'Cell Cycle: G2/M DNA 
Damage Checkpoint Regulation! 'Mitochondrial Dysfunc- 
tion', and 'Role of BRCA2 in DNA Damage Response'. 
The topmost activated function was cell survival' (p-value 
5.08E-03; 89 molecules; activation z-score 5.66). The top- 
most inhibited function was 'organismal death' (p-value 
1.96E-03; 150 molecules; activation z-score -11.38). Similar 
results were obtained using the more limited set of K/P- 
DEGs with the cutoff of >2-fold change. None of the func- 
tion enrichments were related to inflammation or response 
to infection. One possible explanation for this result is 
greater numbers of leukocytes had infiltrated P.a. -infected 
wounds, importing more inflammation transcripts than in 
-infected wounds. These data suggest that after 24 hours 
of active infection, cutaneous cells of K.p. -infected wounds 
activated pathways for surviving the particular conditions of 
the infection such as hypoxia and reactive oxygen species. 

POD4: Cutaneous cells are suggested to activate an 
"Integrated stress response" 

Since the response to infection involves cells infiltrating 
the wound that only contribute transcripts to the wound, 
Dnr-transcripts come from cutaneous cells. The 27%, or 
133, I/U-DEGs that were downregulated to a greater ex- 
tent in K.p.- than P. a. -infected wounds had as their top- 
most enrichment in the category 'Physiological System 
Development and Function' the function annotation 'Hair 
and Skin Development and Function' (p-value 6.20E-06 
to 1.95E-02; 7 molecules), additionally suggesting these 
transcripts come from cutaneous cells. These 133 I/U- 
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DEGs' topmost significant canonical pathway enrichments 
were 'EIF2 Signaling! 'Regulation of eIF4 p70S6K Signaling' 
and mTOR Signaling'. Additionally, the topmost activated 
upstream regulator inferred for the 133 genes was the group 
N-cor (Activation z-score = 2.0), consisting of NCOR1 and 
NCOR2, nuclear receptor co-repressors that mediate basal 
transcriptional silencing [43]. In sum, these results suggest 
that cutaneous cells of K.p. -infected wounds activated an 
"integrated stress response" which involves downregulation 
of general translation and upregulation of specific transcripts 
to survive the stressors of K.p. infection [44,45]. 

The integrated stress response was also suggested for 
P.a. -infected wounds as the top three canonical pathway 
enrichments for 212 unique Dnr DEGs of P. ^.-infec- 
ted wounds were 'Mitochondrial Dysfunction; 'Oxidative 
Phosphorylation' and EIF2 signaling. These as well as 
'Regulation of eIF4 and p70S6K Signaling' and 'mTOR Sig- 
naling' were also the top 5 canonical pathways into which 
enriched 85 Dnr-I/U-DEGs common to K.p.- and P.a. -in- 
fected wounds. Although these data indicate that the cuta- 
neous cells of both types of infected wounds had activated 
an "integrated stress response" to survive the inflamma- 
tory milieu, the P. ^.-unique I/U-DEGs enriched more 
into the 'Cell Death' function category. Using 212 unique 
Dnr-I/U-DEGs to each wound type, 17 from K.p. -infected 
wounds and 72 from P.a. -infected wounds enriched in 
the 'Cell Death' function annotation category. 

POD4: Discordantly regulated transcripts between 
K.p.- and Ra.-infected wounds. 

Of the 9 (1%) discordantly regulated I/U-DEGs between 
the two types of infected wounds, four were significantly 
Upr in P.a. and Dnr in K.p.: SPTAN1, ATP6V1F, 
ATG4A, and EEA1 (K/P-DEGs were 3.8-, 3-, 1.8- and 
1.8-fold higher in P.a., respectively). Each of these genes 
has a function related to intracellular organelles [41]. 
SPTAN1 encodes the filamentous cytoskeletal protein 
spectrin (alpha subunit, aka Fodrin), which functions 
as a membrane-stabilizing scaffold required for the 
function of intracellular organelles. ATP6V1F encodes a 
component of vacuolar ATPase (V-ATPase), which me- 
diates acidification of eukaryotic intracellular organelles. 
ATG4A encodes a cysteine protease required for autoph- 
agy and is implicated in the processing of mature autopha- 
gosomes to facilitate their delivery to the lysosomal 
network. The EEA1 -encoded protein binds phospholipid 
vesicles containing phosphatidylinositol 3-phosphate and 
participates in endosomal trafficking. Notably, P.a. viru- 
lence factor ExoS (expressed by PA01 in our in vitro ex- 
periments, unpublished observation) and in particular its 
ADP-ribosyltransferase activity, is required for intravacuo- 
lar replication of P.a. in epithelial cells [46]. 

The 5 I/U-DEGs that were elevated in K.p.- and sup- 
pressed in P.a. -infected wounds were LMOl, BRMS1L, 



CHD1L (K/P-DEGs were 2.5-, 1.4, and 1.5-fold higher in 
K.p.), as well as two transcribed loci EB374055 and 
DN885016. The three known genes encode chromatin 
binding proteins [41]. The LMOl -encoded protein is 
suggested to regulate transcription by competitively 
binding to specific DNA-binding transcription factors. 
BRMS1L encodes a component of histone deacetylase 
complexes. CHD1L encodes a DNA helicase that plays a 
role in chromatin-remodeling following DNA damage. 
One possibility is that K.p. -infected wounds retained 
capacity to remodel chromatin in response to the stress 
of the infectious inflammation. 

POD4: Network with upstream regulation by type I 
interferons unique to Ra.-infected wounds 

Seeking to identify differences in upstream regulators as- 
sociated with the gene expression patterns of the infected 
wound types, we analyzed upstream activating regulators 
of bacteria. They were markedly similar between the two 
wound types. Given the extensive overlap of identical and 
similar genes, to further search for discriminating features 
of the wound-response to P.a. and K.p. infections, we ana- 
lyzed I/U-DEGs grouped by predicted activated upstream 
regulators (z-score > 2) in the categories, 'Cytokine; 'G -pro- 
tein-coupled receptor; 'Growth factor; 'Ligand-dependent 
nuclear receptor; 'Transmembrane receptor; and 'Other'. 
The I/U-DEGs of each wound type without direct regula- 
tory connections to other molecules in the same group, 
based on the IPA database, were excluded from the ana- 
lysis. There were 89 upstream regulators of the 134 mole- 
cules common to both infected wound types, 69 upstream 
regulators of the 113 molecules unique to P.a. -infected 
wounds, and 3 predicted upstream regulators of the 66 
molecules unique to /^.-infected wounds. 

Between the K.p. -P.a. -common and P.a. -unique I/U-DEGs, 
44 upstream regulators overlapped. The remaining upstream 
regulators summed to 25 and are inferred to regulate 56 
downstream I/U-DEGs unique to P.a. -infected wounds. Fea- 
tures of this network include regulation by type I and type II 
interferons, and 16 downstream molecular targets with more 
than four regulatory connections: INFG, STAT1, STAT4, 
PIM1, PIM2, CCL5, CSF2, MCL1, MMP1, MMP3, IRF1, 
IRF7, RSAD2, EGR1, DHX58, and OASL (Additional file 3). 
Although predicted to be activated by type I interferons, nei- 
ther the type I interferons (IFNA, IFNA2, IFNA5, IFNL3, 
and IFNA14) nor the type I interferon receptor, IFNAR1, 
that were on our microarray were upregulated. 

POD6 wounds have fewer inflammation and stress-response 
transcripts, but Upr-transcripts from Ra.-infected wounds 
form a network with predicted upstream regulation by type 
I interferons as on POD4 

By POD6, wounds infected for 24 hours and then 
treated with antibiotic for two days were likely colonized 
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with biofilm, as approximately 10 5 and 10 4 CFU of K.p. 
and P. a., respectively, were recovered from the wounds 
(our previous study [12]), which is 1 to 2 logs fewer than 
the 10 6 CFU inoculated on POD3. Although about 10- 
fold fewer P.a. CFUs were recovered per wound, P. a. is 
known to replicate in vacuoles of epithelial cells [46]. 
Consistent with low bacterial burden as biofilm, both in- 
fected wound types had far fewer I/U-DEGs on POD6 
than on POD4 (Table 1). 

To assess similarities and differences in the response of 
the wounds to colonization with K.p. vs. P. a. biofilm, we 
analyzed enrichments of the common and the unique I/U- 
DEGs of K.p.- and P. <2. -infected wounds. We analyzed the 
top fold-change I/U-DEGs, 100 Upr and 100 Dnr that 
were significant (adjusted ^-value < 0.05), in order to com- 
pare nearly equal numbers of I/U-DEGs, since I/U-DEGs 
with >2 fold-change were so few. Neither K.p. -P.a. -com- 
mon-I/U-DEGs nor the K.p. -unique I/U-DEGs enriched 
in functions with z-score >2.0. In contrast, the 94 
I/U-DEGs unique to P. <2. -infected wounds enriched in 
functions with predicted activation (z-score 2.7 to 2.0): ac- 
tivation of T lymphocytes; antiviral response; activation of 
cells; activation of antigen presenting cells; cytotoxicity of 
lymphocytes; 'killing of cells; cell death of immune cells; 
activation of leukocytes; activation of cytotoxic T cells; 
antiviral response of cells; response of mononuclear 
leukocytes; 'fragmentation of DNA; cytotoxicity of T lym- 
phocytes; quantity of cytotoxic T cells; maturation of 
phagocytes; maturation of antigen presenting cells; Im- 
mune response of cells; and activation of macrophages'. 
While the 147 I/U-DEGs of K.p. -infected wounds barely 
made any regulatory connection amongst themselves, the 
94 I/U-DEGs of P.a. -infected wounds connected into a 
25-molecule regulatory network in which the major hubs 
were STAT1 and IRF7, with each having 16 direct regula- 
tory connections. This 25-molecule network had several 
genes in common with the upstream-interferon-predicted 
network of I/U-DEGs unique to P. <2. -infected wounds on 
POD4. These common genes between the POD4 and 
POD6 networks are STAT1, IRF7, and DHX58. Addition- 
ally, the POD6 network includes PSMB10, GBP1, and 
OAS2, while the POD4 network includes their highly re- 
lated molecules PSMB9, GBP4, and OASL. These 25 mol- 
ecules' topmost five enrichments into canonical pathways 
are shown in Table 2. The molecule names, and the log ra- 
tios and j?-value of their I/U-DEGs are listed in Additional 
file 4. Furthermore, the Upr I/U-DEGs enriched as their 
topmost network: 'Infectious Disease, Inflammatory Re- 
sponse, and Antimicrobial Response (/?-value lxl0~ 47 right 
tailed Fisher's exact test). Table 3 lists the topmost five up- 
stream mammalian receptors inferred to activate these 
I/U-DEGs. Although type I interferons are prominent 
among these, neither they nor their receptor were upregu- 
lated in our experiments. 



Analysis of the Dnr I/U-DEGs for associations with 
upstream regulators and downstream functions found 
that despite the presence of upregulated inflammatory/ 
antimicrobial transcripts, activation of stress-response 
pathways was not evident as on POD4. 

POD12: More inflammation and inhibited healing 
suggested for Ra.-infected wounds 

On POD 12, uninfected wounds had completely epithelia- 
lized, K.p. -infected wounds were nearly completely epithe- 
lialized (epithelial gap ~1.7 mm), and P.a.- infected 
wounds were more open (epithelial gap ~4 mm) [12]. In 
both K.p.- and P.a. -infected wounds, biofilm was seen on 
POD 12 using scanning electron microscopy [12]. 

The I/U-DEGs common to K.p.- and P.a. -infected 
wounds (Upr and Dnr) enriched in the functions 'Vasculo- 
genesis' and Angiogenesis; predicted to be activated with 
z-scores 2.2 and 2.0 (p-value 1.22E-02 and 1.97E-03, re- 
spectively). The top-scoring function was 'Connective Tis- 
sue Disorders' (p- values 1.27E-06 to 2.33E-03), which 
included collagen type XI, alpha 1J collagen, type XII, 
alpha I', collagen, type V, alpha % 'fibrillin V, 'interleukin 18 
binding protein; 'lysyl oxidase-like % matrix metallopepti- 
dase 1 (interstitial collagenase); 'thrombospondin 2'. All 
were downregulated except IL18BP and MMP1, which 
were upregulated. 

We also analyzed K.p.- versus P. <z. -infected wound dif- 
ferential expression (K/P). Transcripts significantly ele- 
vated in P.a.- relative to /^.-infected wounds summed 
to 118. These molecules significantly enriched into In- 
flammatory Response', 'Connective Tissue Disorders', 'In- 
flammatory Disease', 'Cellular Movement', 'Hematological 
System Development and Function' and 'Immune Cell 
Trafficking' (movement and migration of phagocytes) 
as topmost functions. These 118 genes are shown in 
Additional file 5. Notably, the gene with the largest fold- 
change was PLK5, which was 10-fold elevated in P.a.- as 
compared to K.p. -infected wounds. PLK5 was upregu- 
lated in serum-starved fibroblasts as they exited the 
cell cycle, and PLK5 overexpression arrested cells at a 
GO/Gl-like stage [47]. The topmost predicted upstream 
regulator for the 118 molecules elevated in P. <z. -infected 
wounds was FGF1 (fibroblast growth factor 1, acidic), 
which trended toward an inhibited activation state 
(z-score -0.818; ^-value 3.85E-04). 

Conversely, molecules elevated in K.p.- relative to .Re- 
infected wounds summed to 86 (Additional file 6). A 
single activating upstream regulator was the synthetic 
glucocorticoid methylprednisone (z-score 2.45, ^-value 
8.88E-03, 7 genes: SCAMPI, LDLR, HSD17B12, SDC2, 
CSNK2B, ATIC and FAH). Notably, endogenous gluco- 
corticoid has been suggested to suppress granulation tis- 
sue formation, as mice with a DNA-binding-defective 
glucocorticoid receptor had early wounds with enlarged 
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Table 2 Top five canonical pathway enrichments of 25 molecules from the network on POD6 unique to P.a.-infected 
wounds 



Name 


p-value 


Ratio 


Activation of IRF by Cytosolic Pattern Recognition Receptors 


1.82E-08 


5/72 (0.069) 


Crosstalk between Dendritic Cells and Natural Killer Cells 


1.11E-07 


5/96 (0.052) 


Role of Pattern Recognition Receptors in Recognition of Bacteria and Viruses 


1.54E-07 


5/106 (0.047) 


Role of RIG 1 -like Receptors in Antiviral Innate Immunity 


3.22E-07 


4/49 (0.082) 


Communication between Innate and Adaptive Immune Cells 


6.08E-06 


4/109 (0.037) 



granulation tissue and a high fibroblast density [48]. 
PDGF BB was also a predicted upstream regulator, 
though it did not reach the activating z-score cutoff of 
2.0 (z-score 1.41, Rvalue 1.22E-02, 5 genes: SLC7A7, 
MGST2, THBD, FMOl and LDLR). Additionally, EVPL 
(Envoplakin) was 2.8-fold upregulated in K.p. -infected 
wounds over uninfected wounds and was 2-fold more 
highly expressed in K.p.- than P.a. -infected wounds. 
Envoplakin is a component of the keratinocyte cornified 
envelope, and is suggested to link it to desmosomes and 
intermediate filaments [39]. 

q-RT-PCR validation of microarray results 

Analysis of gene expression using q-RT-PCR of the same 
RNA samples as used for microarray analysis validated 
the microarray results (Additional file 7: Table SI). The 
selected genes represented those that were highly regu- 
lated either in K.p. or P.a. infected wounds. The qRT- 
PCR results confirmed the trend in expression of these 
selected genes revealed by microarray. 

Discussion 

Wound gene expression patterns differed greatly be- 
tween PODs 4, 6 and 12. This was expected as the 
wound states are different on these days, and respect- 
ively involve acute bacterial-induced inflammation, 
stabilization of biofilm, and healing. However, the wound 
responses to K.p. and P.a. infections showed similarities, 
especially on POD4. But there were significant differ- 
ences. On POD4, the Upr gene expression patterns 



indicating inflammation were largely similar between the 
two infected wound types. But P. <2. -infection induced a 
more intense response, with Dnr-transcripts indicating a 
stronger stress-response and the death of cutaneous 
cells. The greater inflammatory response to P.a. infec- 
tion is likely associated with its pathogenic virulence that 
could have resulted in a larger bacterial load as com- 
pared to K.p. infection ([13] and unpublished observa- 
tion). On POD6, after antibiotic killed the majority of 
planktonic bacteria and selected for biofilm, P. <2. -infected 
wounds showed a network of inflammation/infection-re- 
sponse molecules predicted to be regulated by type I in- 
terferons, which was not associated with enrichments 
indicative of a stress response of cutaneous cells. The 
network of inflammation/infection-response molecules 
on POD6 was similar to a network unique to P.a. -in- 
fected wounds on POD4 that was minor part of the ex- 
tensive inflammation/infection-response on that day. On 
POD 12, P. <2. -infected wounds showed indications of 
greater inflammation and impaired healing, which may 
be partly related to a larger P.a. biofilm load at this time 
as compared to K.p. biofilm, as suggested by our previ- 
ous studies [13,49]. 

POD4 

Inflammation/infection-response 

Seeking potential mechanisms of the very similar but 
much more extensive inflammatory host-response of 
P.a. vs. K.p. infection, we analyzed predicted upstream 
activating molecules. Those from bacteria were similar 



Table 3 Topmost five upstream regulators predicted for the 55 Upr l/U-DEGs unique to P.a.-infected wounds on 
POD6 in the categories Cytokine, G-protein-coupled receptor, Growth factor, Ligand-dependent nuclear receptor, 
and Transmembrane receptor 



Upstream 
regulator 


Log 
ratio 


Activation 
z-score 


p-value of 
overlap 


Target molecules 
in dataset 


TLR3 


0.340 


3.111 


7.13E-19 


CXCL10, DDX58, DHX58, GBP4, GCH1, IFIH1, IFIT1, IFIT3, 
IL15, IL23A, IRF7, STAT1 , TLR3, TNFSF10, USP18, ZNFX1 


IFN Beta 




3.395 


1.34E-16 


CXCL10, IFIH1, IFIT1, IFIT1B, IL15, IRF7, STAT1, TAPBP, 
TLR3, TNFSF10, USP18, XAF1 


Interferon alpha 




3.913 


1.67E-16 


CXCL10, DDX58, DHX58, IFIH1, IFIT1, IFIT1B, IFIT3, IL15, 
IRF7, PRF1, STAT1, TAPBP, TLR3, TNFSF10, USP18, WARS 


IFN type 1 




2.087 


3.44E-15 


CXCL10, DDX58, DHX58, IFIH1, IFIT1, IFIT1B, STAT1, TNFSF10, UBA7 


Ifn 




2.793 


4.54E-12 


CXCL10, DDX58, DHX58, IFIH1, IFIT1, IL15, IRF7, STAT1 , TLR3 
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for both infected wound types. To further search for 
specific gene expression patterns that discriminate the 
two infected wound types, we analyzed other upstream 
regulators in categories of cytokines, growth factors and 
signaling receptors. Twenty-five upstream regulators, 
predominated by type I interferons, were inferred to regu- 
late 56 I/U-DEGs unique to P.a. -infected wounds that 
enriched into the functions Immune Cell Trafficking' 
(36 molecules) and Infection (39 molecules; ^-values 
1.24E-26 - 1.04E-7). In the Infection category twice as 
many, 34, enriched into the subcategory 'Viral Infection' 
(p-value 2.56E-18) as enriched into 'Bacterial Infection' 
(17 molecules, p- value 8.37E-13). This network could 
contribute to upregulated IKBKE, which was upregulated 
10-fold in P.a. -infected wounds but only 2-fold in K.p. -in- 
fected wounds. IKBKE is a non-canonical IkappB kinase 
with an important role in regulating antiviral signaling 
pathways. It is activated downstream of the recognition of 
LPS and nucleic acids by their PPRs— TLR3 (dsRNA), 
TLR4 (LPS), IFIH1 (long dsRNA), DDX58 (5'ppp ssRNA, 
short dsRNA), and ZBP1 (dsDNA) [50] — and activated 
IKBKE phosphorylates p65/RelA, interferon response fac- 
tors 3 and 7 (IRF3 and IRF7), and STAT1. Consistent with 
a role in activating IKBKE, the network unique to P.a. -in- 
fected wounds on POD4 includes upregulated TLR3, 
IFIH1, DDX58, IRF3 and STAT1. The network unique to 
P.a. -infected wounds on POD4 and the similar network 
on POD6 could possibly be related to the intracellular rep- 
lication of P.a. [46]. 

Some gene expression induced by P. ^.-infection as 
compared to K.p. -infection probably involves known 
virulence mechanisms of P.a. [46,51], including the type 
III secretion effectors ExoS and ExoT that are expressed 
by the P.a. strain PAOl that we used for our studies 
[52]. A previous study identified genes with induced ex- 
pression dependent on ExoT activity [51] that are some 
of the same molecules that were highly expressed in 
P.a.- as compared to K.p. -infected wounds (e.g., ADM, 
DUSP5 and NR4A1). Additionally, many molecules tar- 
geted by ExoS activities (Ras, Ral, Rabs, and ERM pro- 
teins; as referenced in [53]) were differentially expressed 
in P.a.- vs. K.p. -infected wounds. For example, the ERM- 
like molecule NF2 was one of the 36 genes more upreg- 
ulated in K.p.- than P. <z. -infected wounds. 

Stress-response 

In parallel to the greater number of Upr-transcripts and 
their higher upregulation in P. <2. -infected wounds on 
POD4, P. <2. -infected wounds downregulated more tran- 
scripts, and most of the common Dnr-transcripts were 
downregulated to lower levels as compared to K.p. -in- 
fected wounds. Dnr-transcripts come from cutaneous 
cells, since cells that infiltrate the wound in response to in- 
fection only contribute transcripts. The Dnr-transcripts of 



infected wounds enriched in stress-response pathways 
and networks; these were several-fold more extensive in 
P.a. -infected wounds. 

Given that a constant quantity of RNA was hybridized 
to each microarray, more leukocyte infiltration in P.a.- 
infected wounds would result in more leukocyte RNA 
and proportionally less fibroblast and keratinocyte RNA. 
But this does not appear to fully explain the greater 
extent of Dnr- transcripts in P.a. wounds on POD4 
because P. <2. -infected wounds had 4-fold more Upr- 
transcripts and twice as many, 8-fold, more Dnr-transcripts 
than K.p. -infected wounds (Table 1). And 26% of overlap- 
ping Dnr I/U-DEGs were more downregulated in K.p.- than 
P.a. -infected wounds. Also, stress responses are known to 
involve transcriptional and posttranscriptional processes in- 
cluding the specific stabilization and destabilization of 
mRNA as well as selective translation and transcription, 
and depend on the particular type of stress [54,55]. 

The downregulated transcripts assembled into net- 
works with highly connected nodes such as CAND1 and 
cullins, which are involved in protein degradation, and 
the COP9 signalosome subunits. The COP9 signalosome 
functions in the regulation DNA repair, checkpoint con- 
trol and chromatin remodeling [41]. DNA damage was a 
top activated function enrichment of Dnr I/U-DEGs 
unique to P. <2. -infected wounds, which is a likely conse- 
quence of reactive oxygen species from the strong im- 
mune response to P.a. infection. The topmost predicted 
upstream regulator was MYC. MYC acts as a transcrip- 
tional amplifier, inducing expression of all genes to 
which it is bound [56]. The predicted inhibited state 
of MYC in P. a. -infected wounds is consistent with gen- 
eral suppression of transcription. Also upstream, PPAR 
(peroxisome proliferator-activated receptor) signaling 
was predicted to be in an inhibited state. PPAR signaling 
regulates many proteins in epidermis, which is very ac- 
tive in lipid metabolism, and where all PPAR and LXR 
(liver X receptor) isoforms are expressed, activating 
keratinocyte differentiation and the formation of the per- 
meability barrier [57]. The mTOR pathway was also pre- 
dicted to be in an inhibited state. mTOR is a conserved 
checkpoint protein that senses and integrates growth 
factor signals, and levels of nutrients, oxygen and energy. 
Major downstream targets of mTOR signaling are com- 
ponents of the translation machinery [58]. In contrast, 
RARgamma signaling was predicted to be activated, and 
RXRalpha-RARgamma heterodimers are necessary for 
growth arrest [59]. The Dnr-transcripts of P.a. -infected 
wounds enriched 4-fold more into the cell death func- 
tion annotation than did an equal number of Dnr- 
transcripts of K.p. -infected wounds. 

Previous studies have suggested that eIF2 signaling is 
activated during bacterial infections in vitro [60-63]. In 
epithelial cells, intracellular bacteria have been observed 
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to up-regulate several proinflammatory cytokine genes 
(TNFA, IL1B and IL6) and chemokine genes (IL8 and 
CCL20) (as referenced in [63] and therein). These genes 
were among the most highly expressed genes in our 
P. <z. -infected wounds. Also, human vaginal epithelial 
cells (HVECs) exposed to Staphylococcus aureus strongly 
upregulated chemokines and cytokines [64]. Among 
these, on our microarray, CCL20 was significantly upreg- 
ulated only in P. <2. -infected wounds (5-fold) and CXCL2 
was upregulated 3- and 14-fold, respectively, in K.p.- and 
P. <z. -infected wounds. However, our data cannot discrim- 
inate the source of these Upr-transcripts, leukocytes vs. 
cutaneous cells. 

EIF2 signaling was among the top three canonical 
pathway enrichments for Dnr I/U-DEGs of both infected 
wounds. EIF2 is phosphorylated by several kinases in re- 
sponse to protein disequilibrium in the endoplasmic 
reticulum (ER) and cytoplasm, caused by diverse cellular 
stresses such as nutrient deprivation, exposure to the 
mTOR inhibitor rapamycin, high salt, redox perturb- 
ation, and/or DNA damage [65,66]. EIF2 phosphoryl- 
ation shuts off global translation but allows synthesis of 
specific transcription factors and their target genes that 
support recovery from the stress (e.g., ER chaperones, 
oxidoreductases, and ERAD components). The central 
transcription factors in EIF2 signaling, ATF4 and ATF5, 
were not on our microarray, but the related transcrip- 
tion factors, ATF7IP, ATF6, ATF1, BATF2, ATF3, ATF2, 
and ATF7IP2 were upregulated in P. a.- relative to K.p.- 
infected wounds. 

Beyond a threshold of stress, apoptosis can result. An 
apoptotic switch involving the Unfolded Protein Re- 
sponse (UPR) has been described [67]. The switch to the 
"terminal UPR" is thought to involve transcription of 
the pro-apoptotic transcription factor CHOP (DDIT3), 
which was upregulated 4.6-fold in K.p. -infected wounds 
and 13-fold in P. <2. -infected wounds relative to unin- 
fected wounds. CHOP has been suggested to induce 
death by promoting protein synthesis through increased 
transcription of GADD34 (PPP1R15A), which dephos- 
phorylates EIF2 to reverse the inhibition of translation, 
exacerbating the preexisting ER stress [68,69]. Notably, 
in our study, GADD34 (PPP1R15A) was 30-fold upregu- 
lated in P. <2. -infected wounds, relative to uninfected 
wounds, and was 7.5-fold more highly expressed in 
P. a.- than K.p. -infected wounds. 

POD6 

On POD6, the wounds containing biofilm expressed 
magnitudes fewer inflammation transcripts than before 
antibiotic treatment. Although, 10-fold fewer P. a. than 
K.p. CFU were recovered from wounds on POD6 [12], 
only P. <2. -infected wounds contained Upr I/U-DEGs 
that enriched in inflammation/infection functions. This 



network was virtually the only network found on POD6. 
It was similar to a network on POD4 consisting of Upr 
I/U-DEGs unique to P. <2. -infected wounds, which was a 
minor part of the extensive inflammation/infection-re- 
sponse on POD4. In common between the two networks 
were three identical molecules and three highly similar 
molecules, and both networks contained molecules in- 
volved in nucleic acid sensing with inferred regulation 
by type I interferons. However, type I interferons were 
not upregulated on our microarrays. One possible ex- 
planation for these networks is that they are a conse- 
quence of growth of P. a. within epithelial cells [46] and 
could be regulated via an interferon-independent mech- 
anism such as found in fibroblasts [70,71]. Another 
possibility is that these transcripts could come from a 
specific subpopulation of infiltrating cells such as den- 
dritic cells [72]. Consistent with the greatly reduced in- 
flammation on POD6, the Dnr I/U-DEGs of cutaneous 
cells did not significantly enrich into annotated func- 
tions suggesting a stress response. 

Following two days of antibiotic treatment, POD6 wounds 
differentially expressed fewer genes than on POD4 or 
POD 12. However, a large fraction of the few I/U-DEGs 
present in P. <2. -infected wounds were related to each other 
through regulatory connections, including inferred upstream 
regulation by type I interferons. Such a network of tran- 
scripts could possibly be an identifying marker of P.a. bio- 
film colonization that has potential to develop into chronic 
infection that impairs healing. 

POD12 

Relative to uninfected wounds, both types of infected 
wounds on POD 12 differentially expressed genes that 
enriched into the functions 'Vasculogenesis' and Angio- 
genesis', with predicted activated states, and 'Connective 
Tissue Disorders'. These similar functions between the 
two types of infected wounds are likely related to both 
being open as compared to the uninfected wounds that 
were closed on POD12. 

Transcripts that were more highly expressed in P.a.- 
as compared to K.p. -infected wounds enriched in func- 
tions and networks related to the inflammatory response 
and connective tissue disorders, suggesting that P.a. -in- 
fected wound were more inflamed and disordered. 
Additionally, FGF1 (fibroblast growth factor 1) was the 
topmost upstream regulator of these molecules more 
highly expressed in P. <2. -infected wounds. FGF family 
members possess broad mitogenic and cell survival ac- 
tivities, and FGF1 is a known modifier of endothelial cell 
migration and proliferation, as well as an angiogenic 
factor. However, as an upstream regulator of these 
DEGs, FGF1 was suggested to be in an inhibited state. 
Additionally, PLK5 (polo-like kinase 5) was upregulated 
10-fold in P.a.- relative to -infected wounds. PLK5 
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has been shown to induced Gl arrest when ove rex- 
pressed, to be downregulated in proliferating cells, and 
to accumulate in quiescent serum-starved fibroblast 
[73]. These results could indicate an impaired state of 
healing of P.a. -infected wounds. Conversely, upstream 
regulators of the more highly expressed transcripts in 
K.p.- vs. P.a. -infected wounds included a glucocorticoid 
and PDGF. These predicted upstream regulators may in- 
dicate lower inflammation (glucocorticoid) and active 
healing (PDGF) [74]. Overall these expression data and 
their regulator and function enrichments suggest that 
P. <z. -infected wounds had greater inflammation and K.p.- 
infected wounds had more developed and active healing. 

Conclusions 

These data contribute to defining the immune response of 
wound and proximal skin to bacterial infections of different 
virulence. In the clinical setting, diagnosis of infected 
wounds can be complicated by an absence of signs of in- 
flammation (e.g., as on POD6) and late diagnosis can result 
in impaired healing [75,76]. Current diagnosis of infection is 
based on recoverable CFU which captures only culturable 
bacteria. PCR identifies more bacteria in wound samples, 
but mixed-species infections could have unpredictable viru- 
lence. The host-response counterbalances infections, and 
molecular markers of the host-response to infection may ul- 
timately be useful for diagnosis and treatment. The current 
study suggests that patterns and extents of expression of 
Dnr-transcripts, expressed in cells resident to the wound 
and proximal skin, as well as specific Upr-transcripts, may 
contribute to the prediction or identification of poorly 
healing wounds, and may be useful for indicating wound- 
specific treatments to improve wound outcome. 

Ethics statement 

The animal experiments were performed under animal 
protocol number 2010-1241 approved by Northwestern 
Medical School IACUC. 

Additional files 



Competing interest 

KPL and MJ are employees of the U.S. Government. The work presented is part 
of their official duties. Title 17 U.S.C. §105 provides that 'Copyright protection 
under this title is not available for any work of the United States Government' 
Title 1 7 U.S.C. §101 defined U.S. Government work as work by a military service 
member or employee of the U.S. Government as part of that person's official 
duties. The opinions or assertions contained herein are the private views of 
these authors and are not to be construed as official or as reflecting the views 
of the Department of the Army or the Department of Defense. 

Authors' contributions 

Conceived and designed the experiments: KPL, TC and TAM. Performed the 
experiments: AKS, MRG, WX, SJH, and RDG. Analyzed the data: PD, KPL, TC, 
and WX. Wrote the paper: PD, KPL, and TC. Contributed reagents/materials/ 
analysis tools and helped produce the final manuscript: MJ. Conceived of 
the study, obtained funding, participated in its design and coordination: KPL, 
TAM. Read and approved the final manuscript: KPL, PD, AKS, MRG, MJ, WX, 
SJH, RDG, TC, TAM. 

Acknowledgements 

This work was supported by the US Army Medical Research and Materiel 
Command Military Infectious Diseases (RAD I) and Combat Casualty Care 
(RAD II) Research Directorates. 

Author details 

Microbiology Branch, US Army Dental and Trauma Research Detachment, 
Institute of Surgical Research, 3650 Chambers Pass, Building 3610, JBSA Fort 
Sam Houston, TX 78234, USA. 2 Systems and Integrative Biology, US Army 
Center for Environmental Health Research, Fort Detrick, Frederick, MD, USA. 
3 Division of Plastic Surgery, Feinberg School of Medicine, Northwestern 
University, Chicago, IL 6061 1, USA. 4 The Forsyth Dental Institute, 245 First 
Street, Cambridge, MA 02142, USA. 

Received: 5 November 2013 Accepted: 17 April 2014 
Published: 2 May 2014 

References 

1. Cooper L, Johnson C, Burslem F, Martin P: Wound healing and 
inflammation genes revealed by array analysis of 'macrophageless' PU.1 
null mice. Genome Biol 2005, 6:R5. 

2. Singer AJ, Clark RA: Cutaneous wound healing. N Engl J Med 1 999, 
341:738-746. 

3. Marucha PT, Kiecolt-Glaser JK, Favagehi M: Mucosal wound healing is 
impaired by examination stress. Psychosom Med 1998, 60:362-365. 

4. Sroussi HY, Williams RL, Zhang QL, Villines D, Marucha PT: Ala42S100A8 
ameliorates psychological-stress impaired cutaneous wound healing. 
Brain Behov Immun 2009, 23:755-759. 

5. Liechty KW, Kim HB, Adzick NS, Crombleholme TM: Fetal wound repair 
results in scar formation in interleukin-10-deficient mice in a syngeneic 
murine model of scarless fetal wound repair. J Pediatr Surg 2000, 
35:866-872. discussion 872-863. 

6. Peranteau WH, Zhang L, Muvarak N, Badillo AT, Radu A, Zoltick PW, Liechty 
KW: IL-10 overexpression decreases inflammatory mediators and 
promotes regenerative healing in an adult model of scar formation. 

J Invest Dermatol 2008, 1 28:1 852-1 860. 

7. Menke NB, Ward KR, Witten TM, Bonchev DG, Diegelmann RF: Impaired 
wound healing. Clin Dermatol 2007, 25:19-25. 

8. Pierce GF: Inflammation in nonhealing diabetic wounds: the space-time 
continuum does matter. Am J Pathol 2001, 159:399-403. 

9. Terranova A: The effects of diabetes mellitus on wound healing. 
Plast Surg Nurs 1991, 11:20-25. 

10. Fazli M, Bjarnsholt T, Kirketerp-Moller K, Jorgensen A, Andersen CB, Givskov 
M, Tolker-Nielsen T: Quantitative analysis of the cellular inflammatory 
response against biofilm bacteria in chronic wounds. Wound Repair Regen 
2011, 19:387-391. 

1 1 . Gurjala AN, Geringer MR, Seth AK, Hong SJ, Smeltzer MS, Galiano RD, Leung 
KP, Mustoe TA: Development of a novel, highly quantitative in vivo 
model for the study of biofilm-impaired cutaneous wound healing. 
Wound Repair Regen 201 1 , 1 9:400-41 0. 

12. Seth AK, Geringer MR, Galiano RD, Leung KP, Mustoe TA, Hong SJ: 
Quantitative Comparison and Analysis of Species-Specific Biofilm 



Additional file 1: List of gene probes and associated normalized log2 
fold-changes and adjusted p-values. 

Additional file 2: Genes upregulated to a greater extent in /(.p.-infected 
wounds on POD4. 

Additional file 3: Upstream regulators and the 56 downstream 
molecules unique to P.a.-infected wounds they are inferred to regulate 
on POD4. 

Additional file 4: Regulatory network of P.a.-infected wounds on POD6. 

Additional file 5: Transcripts elevated in P.a- relative to /C.p.-infected 
wounds on POD12. 

Additional file 6: 86 transcripts elevated in K.p.- relative to P.a.-infected 
wounds on POD12. 

Additional file 7: Table SI. RT-qPCR validation of microarray gene 
expression. 



Leung et al. BMC Clinical Pathology 2014, 14:20 
http://www.biomedcentral.eom/1 472-6890/1 4/20 



Page 13 of 14 



Virulence Using an In Vivo, Rabbit-Ear Model. J Am Coll Surg 2012, 
215:388-399. 

13. Seth AK, Geringer MR, Gurjala AN, Abercrombie JA, Chen P, You T, Hong SJ, 
Galiano RD, Mustoe TA, Leung KP: Understanding the host inflammatory 
response to wound infection: an in vivo study of Klebsiella pneumoniae 
in a rabbit ear wound model. Wound Repair Regen 2012, 20:214-225. 

14. Davis SC, Ricotti C, Cazzaniga A, Welsh E, Eaglstein WH, Mertz PM: 
Microscopic and physiologic evidence for biofilm-associated wound 
colonization in vivo. Wound Repair Regen 2008, 16:23-29. 

15. Serralta VW, Harrison-Balestra C, Cazzaniga AL, Davis SC, Mertz M: Lifestyles 
of bacteria in wounds: Presence of biofilms? Wounds 2001, 13:29-34. 

16. Bowler PG: Wound pathophysiology, infection and therapeutic options. 
Ann Med 2002, 34:419-427. 

1 7. Murray CK, Roop SA, Hospenthal DR, Dooley DP, Wenner K, Hammock J, 
Taufen N, Gourdine E: Bacteriology of war wounds at the time of injury. 
Mil Med 2006, 171:826-829. 

18. Aronson NE, Sanders JW, Moran KA: In harm's way: infections in deployed 
American military forces. Clin Infect Dis 2006, 43:1045-1051. 

19. Petersen K, Riddle MS, Danko JR, Blazes DL, Hayden R, Tasker SA, Dunne JR: 
Trauma-related infections in battlefield casualties from Iraq. Ann Surg 
2007, 245:803-811. 

20. Johnson EN, Burns TC, Hayda RA, Hospenthal DR, Murray CK: Infectious 
complications of open type III tibial fractures among combat casualties. 
Clin Infect Dis 2007, 45:409-41 5. 

21. Yun HC, Branstetter JG, Murray CK: Osteomyelitis in military personnel 
wounded in Iraq and Afghanistan. J Trauma 2008, 64:S1 63-1 68. discussion 
S168. 

22. Yun HC, Murray CK, Roop SA, Hospenthal DR, Gourdine E, Dooley DP: 
Bacteria recovered from patients admitted to a deployed U.S. military 
hospital in Baghdad, Iraq. Mil Med 2006, 171:821-825. 

23. Costerton JW: A short history of the development of the biofilm concept. 

In Microbial biofilms. Edited by Ghannoum M, OToole GA. Washington, DC: 
ASM Press; 2004:4-19. 

24. Percival SL, Hill KE, Williams DW, Hooper SJ, Thomas DW, Costerton JW: A 
review of the scientific evidence for biofilms in wounds. Wound Repair 
Regen 2012, 20:647-657. 

25. Schaber JA, Triffo WJ, Suh SJ, Oliver JW, Hastert MC, Griswold JA, Auer M, 
Hamood AN, Rumbaugh KP: Pseudomonas aeruginosa forms biofilms in 
acute infection independent of cell-to-cell signaling. Infect Immun 2007, 
75:3715-3721. 

26. Bjarnsholt T, Kirketerp-Moller K, Jensen PO, Madsen KG, Phipps R, Krogfelt K, 
Hoiby N, Givskov M: Why chronic wounds will not heal: a novel hypothesis. 

Wound Repair Regen 2008, 1 6:2-1 0. 

27. James GA, Swogger E, Wolcott R, Pulcini E, Secor P, Sestrich J, Costerton JW, 
Stewart PS: Biofilms in chronic wounds. Wound Repair Regen 2008, 
16:37-44. 

28. Percival S, Bowler P: Biofilms and their potential role in wound healing. 

Wounds 2004, 16:234-240. 

29. Gurtner GC, Werner S, Barrandon Y, Longaker MT: Wound repair and 
regeneration. Nature 2008, 453:314-321. 

30. Gutierrez-Fernandez A, Inada M, Balbin M, Fueyo A, Pitiot AS, Astudillo A, 
Hirose K, Hirata M, Shapiro SD, Noel A, Werb Z, Krane SM, Lopez-Otin C, 
Puente XS: Increased inflammation delays wound healing in mice 
deficient in collagenase-2 (MMP-8). FASEB J 2007, 21:2580-2591. 

31. Hahm G, Glaser JJ, Elster EA: Biomarkers to predict wound healing: the 
future of complex war wound management. Plast Reconstr Surg 201 1, 
127(Suppl 1):21 S-26S. 

32. Van Delden C, Iglewski BH: Cell-to-cell signaling and Pseudomonas 
aeruginosa infections. Emerg Infect Dis 1998, 4:551-560. 

33. Fang CT, Chuang YP, Shun CT, Chang SC, Wang JT: A novel virulence gene 
in Klebsiella pneumoniae strains causing primary liver abscess and 
septic metastatic complications. J Exp Med 2004, 199:697-705. 

34. Nassif X, Fournier JM, Arondel J, Sansonetti PJ: Mucoid phenotype 
of Klebsiella pneumoniae is a plasmid-encoded virulence factor. 
Infect Immun 1989, 57:546-552. 

35. Schroll C, Barken KB, Krogfelt KA, Struve C: Role of type 1 and type 3 
fimbriae in Klebsiella pneumoniae biofilm formation. BMC Microbiol 2010, 
10:179. 

36. Charles CA, Tomic-Canic M, Vincek V, Nassiri M, Stojadinovic O, Eaglstein 
WH, Kirsner RS: A gene signature of nonhealing venous ulcers: potential 
diagnostic markers. J Am Acad Dermatol 2008, 59:758-771 . 



37. Chen L, Yu G, Langefeld CD, Miller DJ, Guy RT, Raghuram J, Yuan X, 
Herrington DM, Wang Y: Comparative analysis of methods for detecting 
interacting loci. BMC Genomics 201 1, 12:344. 

38. Colwell AS, Longaker MT, Peter Lorenz H: Identification of differentially 
regulated genes in fetal wounds during regenerative repair. Wound 
Repair Regen 2008, 16:450-459. 

39. Smyth GK: Limma: linear models for microarray data. In Bioinformatics and 
Computational Biology Solutions using R and Bioconductor. New York: Springer; 
2005:397-420. [Gentleman R, Carey V, S. D, Irizarry R, Huber W (Series Editor). 

40. Kramer A, Green J, Pollard J Jr, Tugendreich S: Causal analysis approaches 
in Ingenuity Pathway Analysis. Bioinformatics 2014, 30:523-530. 

41. Magrane M, Consortium U: UniProt Knowledgebase: a hub of integrated 
protein data. Database (Oxford) 201 1, 201 1:bar009. 

42. Hamidi T, Algul H, Cano CE, Sandi MJ, Molejon Ml, Riemann M, Calvo EL, 
Lomberk G, Dagorn JC, Weih F, Urrutia R, Schmid RM, lovanna JL: Nuclear 
protein 1 promotes pancreatic cancer development and protects cells 
from stress by inhibiting apoptosis. J Clin Invest 2012, 122:2092-2103. 

43. Mottis A, Mouchiroud L, Auwerx J: Emerging roles of the corepressors 
NCoR1 and SMRT in homeostasis. Genes Dev 2013, 27:819-835. 

44. Harding HP, Zhang Y, Zeng H, Novoa I, Lu PD, Calfon M, Sadri N, Yun C, 
Popko B, Paules R, Stojdl DF, Bell JC, Hettmann T, Leiden JM, Ron D: 

An integrated stress response regulates amino acid metabolism and 
resistance to oxidative stress. Molecular cell 2003, 1 1 :61 9-633. 

45. Teske BF, Fusakio ME, Zhou D, Shan J, McClintick JN, Kilberg MS, Wek RC: 
CHOP induces activating transcription factor 5 (ATF5) to trigger 
apoptosis in response to perturbations in protein homeostasis. 

Mol Biol Cell 201 3, 24:2477-2490. 

46. Hritonenko V, Evans DJ, Fleiszig SM: Translocon-independent intracellular 
replication by Pseudomonas aeruginosa requires the ADP-ribosylation 
domain of ExoS. Microbes Infect 2012, 14:1366-1373. 

47. de Career G, Manning G, Malumbres M: From Plk1 to Plk5: functional 
evolution of polo-like kinases. Cell cycle 201 1, 10:2255-2262. 

48. Grose R, Werner S, Kessler D, Tuckermann J, Huggel K, Durka S, Reichardt 
HM, Werner S: A role for endogenous glucocorticoids in wound repair. 
EMBO reports 2002, 3:575-582. 

49. Seth AK, Geringer MR, Gurjala AN, Hong SJ, Galiano RD, Leung KP, Mustoe 
TA: Treatment of Pseudomonas aeruginosa biofilm-infected wounds with 
clinical wound care strategies: a quantitative study using an in vivo 
rabbit ear model. Plast Reconstr Surg 2012, 1 29:262e-274e. 

50. Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M: KEGG for integration 
and interpretation of large-scale molecular data sets. Nucleic Acids Res 
2012, 40:D109-1 14. 

51. McMorran B, Town L, Costelloe E, Palmer J, Engel J, Hume D, Wainwright B: 
Effector ExoU from the type III secretion system is an important 
modulator of gene expression in lung epithelial cells in response to 
Pseudomonas aeruginosa infection. Infect Immun 2003, 71:6035-6044. 

52. Vance RE, Rietsch A, Mekalanos JJ: Role of the type III secreted 
exoenzymes S, T, and Y in systemic spread of Pseudomonas aeruginosa 
PA01 in vivo. Infect Immun 2005, 73:1 706-1 71 3. 

53. Soong G, Parker D, Magargee M, Prince AS: The type III toxins of 
Pseudomonas aeruginosa disrupt epithelial barrier function. J Bacteriol 
2008, 190:2814-2821. 

54. Murray Jl, Whitfield ML, Trinklein ND, Myers RM, Brown PO, Botstein D: 
Diverse and specific gene expression responses to stresses in cultured 
human cells. Mol Biol Cell 2004, 15:2361-2374. 

55. Kawai T, Fan J, Mazan-Mamczarz K, Gorospe M: Global mRNA stabilization 
preferentially linked to translational repression during the endoplasmic 
reticulum stress response. Mol Cell Biol 2004, 24:6773-6787. 

56. Nie Z, Hu G, Wei G, Cui K, Yamane A, Resch W, Wang R, Green DR, 
Tessarollo L, Casellas R, Zhao K, Levens D: c-Myc is a universal amplifier of 
expressed genes in lymphocytes and embryonic stem cells. Cell 2012, 
151:68-79. 

57. Schmuth M, Jiang YJ, Dubrac S, Elias PM, Feingold KR: Thematic review 
series: skin lipids. Peroxisome proliferator-activated receptors and liver 
X receptors in epidermal biology. J Lipid Res 2008, 49:499-509. 

58. Hay N, Sonenberg N: Upstream and downstream of mTOR. Genes Dev 
2004, 18:1926-1945. 

59. Germain P, Chambon P, Eichele G, Evans RM, Lazar MA, Leid M, 
De Lera AR, Lotan R, Mangelsdorf DJ, Gronemeyer H: International 
Union of Pharmacology. LX. Retinoic acid receptors. Pharmacol Rev 
2006, 58:712-725. 



Leung et al. BMC Clinical Pathology 201 4, 14:20 Page 1 4 of 1 4 

http://www.biomedcentral.eom/1 472-6890/1 4/20 



60. Shrestha N, Bahnan W, Wiley DJ, Barber G, Fields KA, Schesser K: Eukaryotic 
initiation factor 2 (elF2) signaling regulates proinflammatory cytokine 
expression and bacterial invasion. J Biol Chem 2012, 287:28738-28744. 

61. Subrahmanyam YV, Yamaga S, Prashar Y, Lee HH, Hoe NP, Kluger Y, Gerstein 
M, Goguen JD, Newburger PE, Weissman SM: RNA expression patterns 
change dramatically in human neutrophils exposed to bacteria. 

Blood 2001,97:2457-2468. 

62. Bohn E, Muller S, Lauber J, Geffers R, Speer N, Spieth C, Krejci J, Manncke B, 
Buer J, Zell A, Autenrieth IB: Gene expression patterns of epithelial cells 
modulated by pathogenicity factors of Yersinia enterocolitica. 

Cell Microbiol 2004, 6:129-141. 

63. Li X, Fusco WG, Seo KS, Bayles KW, Mosley EE, McGuire MA, Bohach GA: 
Epithelial Cell Gene Expression Induced by Intracellular Staphylococcus 
aureus. Int J Microbiol 2009, 2009:753278. 

64. Peterson ML, Ault K, Kremer MJ, Klingelhutz AJ, Davis CC, Squier CA, 
Schlievert PM: The innate immune system is activated by stimulation of 
vaginal epithelial cells with Staphylococcus aureus and toxic shock 
syndrome toxin 1. Infect Immun 2005, 73:2164-2174. 

65. Narasimhan J, Staschke KA, Wek RC: Dimerization is required for activation 
of elF2 kinase Gcn2 in response to diverse environmental stress 
conditions. J Biol Chem 2004, 279:22820-22832. 

66. Peidis P, Papadakis Al, Muaddi H, Richard S, Koromilas AE: Doxorubicin 
bypasses the cytoprotective effects of elF2alpha phosphorylation and 
promotes PKR-mediated cell death. Cell Death Differ 201 1, 18:145-154. 

67. Papa FR: Endoplasmic reticulum stress, pancreatic beta-cell degeneration, 
and diabetes. Cold Spring Harbor perspectives in medicine 2012, 2:a007666. 

68. Baird TD, Wek RC: Eukaryotic initiation factor 2 phosphorylation and 
translational control in metabolism. Adv Nutr 2012, 3:307-321. 

69. Han J, Back SH, Hur J, Lin YH, Gildersleeve R, Shan J, Yuan CL, Krokowski D, 
Wang S, Hatzoglou M, Kilberg MS, Sartor MA, Kaufman RJ: ER-stress-induced 
transcriptional regulation increases protein synthesis leading to cell death. 
Nat Cell Biol 2013, 15:481-490. 

70. Hasan M, Koch J, Rakheja D, Pattnaik AK, Brugarolas J, Dozmorov I, Levine B, 
Wakeland EK, Lee-Kirsch MA, Yan N: Trexl regulates lysosomal biogenesis 
and interferon-independent activation of antiviral genes. Nat Immunol 
2013, 14:61-71. 

71. Dixit E, Boulant S, Zhang Y, Lee AS, Odendall C, Shum B, Hacohen N, 
Chen ZJ, Whelan SP, Fransen M, Nibert ML, Superti-Furga G, Kagan JC: 
Peroxisomes are signaling platforms for antiviral innate immunity. 
Ce//2010, 141:668-681. 

72. Clavarino G, Claudio N, Dalet A, Terawaki S, Couderc T, Chasson L, Ceppi M, 
Schmidt EK, WengerT, Lecuit M, Gatti E, Pierre P: Protein phosphatase 1 
subunit Ppp1r15a/GADD34 regulates cytokine production in 
polyinosinicpolycytidylic acid-stimulated dendritic cells. Proc Natl Acad 
SciUSAlOU, 109:3006-3011. 

73. de Career G, Escobar B, Higuero AM, Garcia L, Anson A, Perez G, Mollejo M, 
Manning G, Melendez B, Abad-Rodriguez J, Malumbres M: Plk5, a polo box 
domain-only protein with specific roles in neuron differentiation and 
glioblastoma suppression. Mol Cell Biol 201 1, 31:1225-1239. 

74. Pierce GF, Tarpley JE, Tseng J, Bready J, Chang D, Kenney WC, Rudolph R, 
Robson MC, Vande Berg J, Reid P: Detection of platelet-derived growth 
factor (PDGF)-AA in actively healing human wounds treated with 
recombinant PDGF-BB and absence of PDGF in chronic nonhealing 
wounds. J Clin Invest 1995, 96:1336-1350. 

75. Nakagami G, Asada M, Sanada H: Innovations in wound infection 
management. Wounds International 2012, 3:4-6. 

76. Asada M, Nakagami G, Minematsu T, Nagase T, Akase T, Huang L, Yoshimura 
K, Sanada H: Novel biomarkers for the detection of wound infection by 
wound fluid RT-PCR in rats. Exp Dermatol 2012, 21:1 18-122. 



doi:1 0.1 1 86/1 472-6890-1 4-20 

Cite this article as: Leung et al.: Dermal wound transcriptomic responses 
to Infection with Pseudomonas aeruginosa versus Klebsiella pneumoniae 
in a rabbit ear wound model. BMC Clinical Pathology 2014 14:20. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



o 



BioMed Central 



